Distributions of inherent structure energies during aging 
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We perform extensive simulations of a binary mixture Lennard- Jones system subjected to a tem- 
perature jump in order to study the time evolution of fluctuations during aging. Analyzing data 
from 1500 different aging realizations, we calculate distributions of inherent structure energies for 
different aging times and contrast them with equilibrium. We find that the distributions initially 
become narrower and then widen as the system equilibrates. For deep quenches, fluctuations in the 
glassy system differ significantly from those observed in equilibrium. Simulation results are par- 
tially captured by theoretical predictions only when the final temperature is higher than the mode 
coupling temperature. 
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When a liquid is rapidly brought out of equilibrium, for 
example by an abrupt change in temperature T, it starts 
to age. The material properties change with the aging 
time t w in the attempt to recover equilibrium at the new 
bath temperature Tb. This aging process can be visual- 
ized as a path in configurational space, that starts from a 
typical equilibrium configuration before the quench and 
which, at long times, should converge to an equilibrium 
configuration characteristic of Tb. When Tb is below the 
glass transition temperature, the aging dynamics never 
stops in the experimentally accessible time window, but 
becomes slower and slower on increasing t w . 

In recent years, several studies have revisited the 
physics of aging 0, @j S 13 ) attempting to associate the 
aging dynamics with a progressive decrease of the fictive 
(or effective) temperature Tf [H, IS IS El • In these stud- 
ies, the configuration of the aging system at time t w is 
associated with a typical configuration explored by the 
equilibrium liquid at Tf. In this case, the aging dynam- 
ics can be modeled as a progressive thermalization of the 
system, quantified by the progressive decrease of Tf . 

A convenient framework for analyzing the thermody- 
namics of supercooled liquids and aging dynamics in nu- 
merical simulations is provided by the inherent structure 
(IS) formalism {js^. In this framework, each configura- 
tion is associated with its closest local minimum, named 
IS, on the potential energy surface (PES) of the liquid. 
A formal thermodynamic description of the system can 
be developed with a proper modeling of the statistical 
properties of the PES, i.e., the distribution of the energy 
minima eis and of their basins of attraction |lfjl ITU . It 
has been shown that in equilibrium supercooled states, 
the probability density P(eis,T) of exploring a mini- 
mum of depth eis is approximately Gaussian, with a 
T-independent variance and with an average eis which 
decreases monotonically with T fl^L liH , fl4| . In the IS ap- 
proach, the aging dynamics can be visualized as the pro- 
gressive exploration of deeper and deeper minima, and a 
connection between aging and equilibrium configurations 
can be made by comparing properties of the explored 



basins. 

Experimental studies of glass forming materials 0, 0] 
and recent numerical simulations |16( have called atten- 
tion to the fact that aging systems, especially for deep 
quenches, may not be associated with a single Tf , i.e., can 
not be uniquely associated with a liquid configuration, 
calling for a further development of the present theories. 
In this Letter, we bring the level of comparison between 
equilibrium and aging to a much more detailed level, by 
studying fluctuations around the average properties. We 
show that the fluctuations frozen in the glass depend on 
the thermal history and are significantly different from 
the fluctuations characteristic of the liquid state. 

We perform molecular dynamics (MD) simulations of 
N = 1000 particles in a cubic box of length L = 9.4, 
for the well characterized Kob- Andersen binary-mixture 
Lennard-Joncs (LJ) model [13 , for which the mode cou- 
pling transition has been estimated at T x = 0.435. We 
employ the Nose-Hoover thermostat with parameters 
chosen so that the kinetic energy thermalizes in a time 
between 300 and 1000 time steps. Time is reported here 
in number of MD time steps, each of which is O.Olio [Tsj - 
We calculate the probability density P(ejs, t w ) of finding 
the liquid at a given time i w in a basin of energy eis > &s 
the system responds to a sudden lowering in T. To this 
end, we perform an ensemble of 1500 independent MD 
runs, each for 10 6 steps. The temperature is switched 
from the initial equilibrium temperature T\ to the new 
T> at t w = 0. For each MD run we select configurations 
at 35 different i w , logarithmically spaced. We perform a 
conjugate-gradient quench on each configuration to ob- 
tain the associated IS. Diagonalization of the Hessian 
matrix evaluated at the IS provides the eigenfrequencies 
{uJi}. Thus for each studied i w we produce an ensem- 
ble of 1500 IS's. We study five different T jumps, with 
sets of (T -> T b ) including (0.55 0.466), (0.6 -> 0.1), 
(0.8 -> 0.25), (0.8 -> 0.446) and (0.8 -► 0.6). 

In Fig. Qwe plot P(ei S ,i w ) for the (0.8 -> 0.446) run, 
as well as the corresponding equilibrium distributions 
P(eis,T) at T and Tb. The figure shows that distribu- 
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FIG. 1: Distributions of eis during aging. Histograms (of unit 
area) showing P(eis,t w ) (filled diamonds) for an ensemble of 
aging systems with T% = 0.8 and T b = 0.446 with t v = 542, 
3405, 7482, 21372, 46954, 174338 and 841500 from right to 
left. Also shown are equilibrium distributions P(eis, T) (open 
circles) for T = 0.8, 0.55 and 0.446 from right to left. 
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tions initially become narrower (equivalently taller), then 
broader again on approaching equilibrium. Fig. [5] shows 
the average e ag and the standard deviation cr ag of the dis- 
tributions as functions of t w . In the cases where T, > T x , 
the variance, which by construction at i w = coincides 
with the equilibrium variance, first decreases and then re- 
laxes towards the equilibrium value. For Tb < T x , within 
the studied time window, we only observe the narrowing. 
We also evaluate the skewness, though we can not discern 
any trend outside of statistical noise. 

Next we attempt to compare the fluctuations observed 
in aging with theoretical predictions. In equilibrium, the 
Helmholtz free energy per particle of the liquid can be 
written as |IcLll^. 



F = eis - TS c (eis) + /vib(eis,T), 



(1) 



where / v ib(eis,T) is the free energy of a basin, av- 
eraged over all basins of depth eis, and S'c(eis) = 
N~ 1 kn hi ^(eis) is the configurational entropy, fi(eis) 
counts the number of basins of energy eis ■ In the 
Gaussian-harmonic approximation (GHA) |l9| . one as- 
sumes that S c and f V ib are quadratic in eis- In the GHA, 
P(eis, T) is Gaussian with a T-independent variance, i.e., 



P(eis,T) 



y/27 



■. exp 



(eis - eis(T)) 2 



2o\ 



(2) 



where the T dependence of eis and the variance a\ are 
expressed in terms of the statistical properties of the PES 
(the variance a 2 and the average energy T of fi(eis)) 
and of the parameters connecting depth and shape of 
the basins (b and c) More precisely, e"is(T) = 

{E - ba 2 )/(l + 2ca 2 c ) - T- 1 a c 2 /(l + 2ccr c 2 ), and cr 2 , = 
N^al/il + 2c<7 2 c ). Below T = 0.6, the GHA provides a 



FIG. 2: Time evolution of the average e ag (a) and of the 
standard deviation a ag (b) of the inherent structure energy 



distribution P(eis,in 
Ti and Tb values. 



Legend indicates the corresponding 



good description of the numerical data [l2l ITl ] , enabling 
us to estimate cr 2 , Eq, b and c. 

In the non-equilibrium case, recent theoretical ap- 
proaches suggest that the free energy can be written as, 



F = eis - TfS'c(eis) + /vib(eis, T h ), 



(3) 



where 5* c is weighted by T f |fj,|a|2J,|2J. The introduction 
of Tf modifies Eq.|21to yield a non-equilibrium probability 
density, in the GHA, 
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27TCT 2 



P(e ls ,T f (t w ),T h ) = 



where e ag and a a g are given by 



: exp 



(eis 



2ct 2 



(4) 



E Tt - al - a 2 c T h b 
2ca 2 T h + Tf ' 



and 



(5) 



(6) 



Note that the theory predicts a gradual increase of the 
variance on increasing i w , i.e., on approaching equilib- 
rium (decreasing Tf). 

The theoretical predictions for the relation between e ag 
and cr a g, parametric in Tf, are shown in Fig. IHIa). In 
Fig-EIb) we show the evolution of the same quantities as 
calculated from the aging runs. In the simulation data, at 
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FIG. 3: Relationship between o" ag and e ag during aging. The 
legend indicates 71 and Tb. Panel (a) shows the theoretical 
predictions based on the GHA, using the landscape parame- 
ters calculated fitting equilibrium data |2l|. Panel (b) shows 
MD data. Filled circles indicate equilibrium values. 



short £ w , the system possesses by construction a distribu- 
tion of IS energies identical to that of the starting T, and 
hence it is not surprising that the numerical results ini- 
tially differ from the value predicted by Eq. El For longer 
times, i.e., for lower eis values, a region where the vari- 
ance increases is indeed observed, but only for quenches 
above T x . The magnitude of change in the variance pre- 
dicted by the theory is significantly smaller than the one 
observed in simulation. Even in the case (0.55 — ► 0.466), 
where the starting T is such that the GHA should hold 
the best, the observed narrowing falls outside the un- 
certainty arising from modeling the equilibrium distribu- 
tions, as well as from numerical errors. 

An implicit assumption in the theory is that, dur- 
ing aging, the liquid explores configurations typical of 
equilibrium |24j and hence that the landscape param- 
eters derived from the equilibrium study can be used 
to calculate the out-of-equilibrium distributions. Previ- 
ous work has shown that this is not necessarily so [25| . 
Here we analyze our extremely large database of inher- 
ent structures in equilibrium and in aging to address 
the question of whether the basin curvature during ag- 
ing differs from the equilibrium case. The quantity 
M = N^ 1 y^f^ -3 In constitutes a good indicator 
of the curvature of a single basin in the harmonic ap- 



FIG. 4: Vibrational properties of basins as a function of 
ds. Panel (a) shows (M) eq (eis) (filled circles, bars in- 
dicating standard deviations) and aging ensemble averages 
(M) (eis). For the dot-dash line we use a binning procedure 
to determine M(eis) from all basins regardless of t v during 
aging for the (0.8 — ► 0.25) run. Panel (b) shows the equilib- 
rium P(M) for T = 0.55 (solid line), and for the case where 
all basins with energy eis(T = 0.55) ± 0.002 are used (short 
dash); and for the aging runs (0.8 — » 0.446) (long dash) and 
(0.8 — » 0.25) (dot dash), also binning the energy. 



proximation [2J, |2£ 



ues (M) eq (ei S ) and (M) ag (e IS 



In Fig. 01a) we plot average val- 
for equilibrium and ag- 
ing runs. This plot shows that for all aging experiments 
with Tb > T x , the eis dependence of M is the same in 
equilibrium and in aging, and hence that the aging liquid 
explores basins the shape of which is the same as the one 
explored in equilibrium. In contrast, for Tb < T x , the 
aging system explores basins of different shape. In the 
(0.8 — > 0.25) case, the system begins to explore basins of 
higher M relative to equilibrium values near eis ~ —7.6. 
In the deep quench (0.6 — ► 0.1) case, the system stays 
in basins of roughly constant / v ib- It is as though the 
system explores basins of lower energy, but not of corre- 
spondingly lower /vib- To discern the origins of the dif- 
ferences in basin sampling in aging and equilibrium, we 
calculate the distribution of M, P(M), for four different 
cases (Fig. EJb)). More specifically, we want to assess if 
differences in (M), at the same average eis value, arise 
from differences in the distributions of sampled eis or 
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if they arise from intrinsic differences in sampled basin 
shapes. For this reason we show, together with the equi- 
librium case for T — 0.55 (for which eig = —7.6228), 
the P{M) for basins of depth e IS = -7.6228 ± 0.002, ex- 
plored in the shallow (0.8 -> 0.446) and deep (0.8 -> 0.25) 
quenches (independent from t w ) and basins explored in 
equilibrium (independent from T). Data in Fig. 0{b) 
show that the equilibrium P{M) coincides with the dis- 
tribution evaluated from all sampled basins with depth 
equal to ejs = —7.6228. It also coincides with the dis- 
tributions evaluated during aging in the Tb > T x runs. 
But it does not coincide with the distributions evaluated 
during aging in the Tb < T x case. Hence, differences be- 
tween equilibrium and aging can not be ascribed to the 
the non- linear dependence of M on eis- We conclude that 
in aging runs with Tb < T x the system explores narrower 
basins (higher M) as compared to equilibrium. 

In conclusion, the results reported in this Letter pro- 
vide high quality data for the evolution of the energy 
distribution in aging, for several different aging scenar- 
ios, and hence provide an important starting poin t for a 
detailed master equation description [2^, l28l E^ ] of 
the dynamics in configuration space. In our simulations, 
when Tb > T x , we find that the distribution initially nar- 
rows during aging, broadening again as the system re- 
equilibrates. We also find that the system samples con- 
figurations which are typically explored in equilibrium. 
In this T range, a description of the aging dynamics via 
a master equation is feasible, provided equilibrium prop- 



erties are fully characterized. By contrast, when T> < T x , 
we only observe the narrowing of the distribution. Fur- 
thermore, the vibrational states explored during aging 
differ from the equilibrium ones. We note that the out-of- 
equilibrium theory, which is built on the hypothesis that 
states explored in aging are the same as those explored 
in equilibrium, predicts only broadening associated with 
the approach to equilibrium. Hence a theory based on 
a single Tf may only apply to the Tb > T x case. We 
also stress that the results reported in this Letter sug- 
gest that fluctuations in arrested states (which on the 
time scale of simulation correspond to Td < T x ) depend 
on the previous thermal history [2^. Glasses with the 
same average values may strongly differ in their frozen 
fluctuations. For deep quenches, the fluctuations frozen 
in the glass are significantly different from the fluctu- 
ations experienced in the liquid state, a feature which 
could explain the cross-over effect HEHEIl and the fail- 
ure of the theories based on one Tf . Finally, we note that 
the decrease of the variance at short times (as shown in 
Fig- El 7 suggests that the initial aging dynamics — which 
may well be the only one accessible on the observation 
time scale for deep quenches — acts in the direction of 
increasing the differences between the equilibrium liquid 
and glass distributions. 
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